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Abstmct-A novel method is used to design wavelets matched 
to the ECG signal. The designed matched wavelet is then used 
for detecting the R peaks of ECG signal through a new detection 
algorithm. No pre-processing or filtering is applied; the signals 
are directly analyzed using the proposed R peak detection 
algorithm. The algorithm is evaluated on real data from MIT- 
BIH Arrhythmia database. Experimental results illustrate that 
the proposed method is faster than the previous methods while 
having similar accuracy. 

Ks^lMinis - Wavelet transform, matched wavelet, ECG signal, 
QRS complex, R peak detection. 

I, INTRODUCTION 

Detection of R peak is the first step in processing an 
ECG signal as it is usually used as a reference within the 
caniiac cyclic and for beat alignment. So far, a vast number 
of algorithms have been proposed for the detection of the 
QRS complex which can be categorized as: 1) syntactic, 2) 
non-syntactic, and 3) hybrid [1], [2]. 

The syntactic algorithms are time consuming and are not 
appnapriate for online applications. In the non-syntactic 
appnaach, the signal is first filtered by a band pass filter or a 
matched filter in order to suppress the T and P waves and the 
noise. The output of the filter is then passed through a 
nonlinear transformer to enhance the QRS complex. But 
different subjects and even different beats of the same subject 
have different QRS frequency bands and also the frequency 
bands of QRS and noise may overlap. 

Wavelet transform is a powerful tool for analyzing non- 
stationary signals, such as ECG. It pnavides a time- frequency 
analysis of the signals by decomposing them into orthogonal 
frequency channels that have the same bandwidth on a 
logarithmic scale. These 'building blocks' are well localized 
both in time and frequency .This feature of the wavelet 
transform is very effective for feature extraction and 
classification [3], [4]. Using the wavelet transform, the 
variations in frequency components of an ECG signal can be 
followed. Thus the accuracy and reliability of the ECG 
feature extraction improves significantly. 

Wavelet-based methods for R peak detection have 
mainly used quadratic spline and Daubechies wavelets which 
have lead to excellent results and robustness against noise 
and baseline drift. Many different wavelet bases with 
different characteristics have been designed and can be used 
for ECG processing [4]-[7]. In addition, it is possible to 
directly design the appnopriate (matched) wavelet for a 
particular application. Matched wavelets can be designed 
such that the error between the original signal and some finite 
resolution wavelet representation of it is minimized [4]-[6]. 

This paper is organized as follows. In Section II, we 
present an introduction to the wavelet transform and then 
explain the design of matched wavelets. In this section, we 



also present a new fast and reliable approach for R peak 
detection which is compatible with different wavelet bases. 
In Section III, through experimental studies, we show the 
power of the new designed matched wavelet and also the new 
approach for R peak detection using the standard, annotated 
MIT-BIH Arrhythmia ECG database. Section IV is our 
concluding remarks. 

II. METHODOLOGY 

The wavelet transform (WT) of a signal f(x) is a 
decomposition of the signal as a combination of a set of basis 
functions [1]. The basis functions are the dilated and 
translated versions of a single prototype wavelet function 

¥(t). So Wj(x), the WT of the signal, is defined as: 

_+« 
Wj(x)=f(x)*lf/Jx)= S-' f f{t)i//{^^^^)dt (1) 

— CO 

The greater the scale factor s, the wider the basis 
function which corresponds to the lower frequency 
components of the signal and vice versa. If we choose the 

scahng factor s such that s = 2^ He Z), then the WT is 
called dyadic WT. Regarding Mallat algorithm, the dyadic 
WT of a discrete-time signal (DWT) is equivalent to an 
octave filter bank and thus can be calculated as follows: 
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S^ 



is a snKothing operator where S^ is equal with the 



digital signal to be analyzed and /2j. , g^. are the coefficients 

of the low pass and high pass FIR filters H(w) and G(w), 

respectively. 5* /(«) can be interpreted as the 

approximation of the original signal at they"' scale whereas 

W f{n) demonstrates the details of the signal at that scale. 

The hi^ pass and low pass filters H (w) and G (w) 
depend on the selected wavelet prototype function Let 

[hQ,h^,...,hj^_^] and [go,gi,..vgw_i ] be the inpjlse 

response of the low pass and high pass filters, respectively. 
The following corxlitions aisure orthogonality of the 
transform so that no information is lost in the deconposition 
process and the filter bank is said to aijqy a perfect 
reconstmction (PR) property [8]. 

Z ^A.2.=L ...... (4a) 
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g„ = {-ir'h 



N-n-l 



,n = ^\,...,N-\ 



(4b) 



2) Designing Matched Wavelet 



The objective function of the methods proposed in [5]- 
[7] is the squared approximation error between the signal of 
interest and its approximation at some finite prescribed scale. 
Since the approximation error is not generally an explicit 
function of the design parameters, an approximation of this 
objective function is required so that design problem is 
converted to a conventional optimization problem [5] - [7]. 

The matching algorithm proposed in [5] is sub-optimal in 
the sense that it is performed on the spectrum magnitude and 
phase independent of one another The algorithm proposed in 
[6] has a more accurate approximation of the squared error 
than the algorithm presented in [7], so this algorithm is 
chosen among the proposed algorithms in [5]-[7] for 
designing wavelets matched to the ECG signal. 

For the given signal /(Itj, the prescribed scale (resolution) 
J, and the dilation factor M, the algorithm seeks for the 
optimal M-band wavelet multi- resolution that represents /fit) 
at scale J. The optimality is measured with respect to 
minimization of frequency norm of the approximation error 
In this work, 2-band wavelet analysis and synthesis is 
implemented (M=2). 

M-band compactly supported orthonormal wavelet bases 
can be uniquely characterized by a unitary sequence known 
as the unitary scaling vector, ho- For the N-length scaling 
vector, the following linear and quadratic constraints exist: 



Af-l 



Y,K{k) = 42 



k=a 



Y,h^{k)\{k + Ml) = 5{l) 



(5) 



(6) 



jt=0 



G iven a length N and the scaling vector ho , there exists a 
unique, compactly supported scaling function 

^n (0 e L^ {R) that satisfies the following lecuision: 



N-\ 



(p,{t) = ^Y.Kik)ni'^t-k) 



(7) 



t-o 



Due to multi-resolution theorem, any function f(t) 
e L^{R) can be miiqudy represented in the following form 

f(t)=^cj^(k)2^^'"^,(y't-k}+^±dj(k),^j,(t) (8) 



^ j=ji} 



Where 



!^,,(?) = 2^'V(2^?-^) 



(9) 

Therefore, all functions f(t) can be practically 
approximated at scale J as 

f{t)^Y.^Ak)2'"ni'^'t-k) (10) 

k 

SO the ^jproximation error can be defined as 
Qfit) = fit)-Y,Cjik)2'"<p,i2't-k) (11) 

Using the parameterization method proposed in [6], the 
scaling vector hg of length N=MK (M=2), can be 
parameterized by K-1 unit vectors v,. Each unit vector itself 



is parameterized by M-1 angle parameters 9. j. , i=0,... , K-1; 

k=0, i, ... , M-2. Following the mathematical steps and 
lemmas in [6], the design of optimal multi-resolution 
analysis, with the Lp error norm, is converted to the following 
optimization problem 



\dw 
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(12) 



wdiere q)^\% the Fourier transform of ^j (0- 

Qice the optimization problem stated in (12) is solved, 
the scaling vector, the scaling function and thus the wavelet 
basis is uniquely obtained 

3) Designing Wavelet Matched to ECG Signal 

In order to preserve compatibility with the previous 
related works, MIT-BIH Arrhythmia, which is a standard 
annotated ECG database, is selected as the data set both for 
training and testing the proposed algorithm. This database 
contains 48 real ECG records which are 30 minutes long and 
have been sampled at 360 Hz. 

A typical cardiac cycle is selected from the first record of 
this database (record#100). To follow the recommendations 
of the AAMI (Association for the Advancement of Medical 
Instrument), the first 5 minute of this record is removed and 
the first cardiac cycle appearing after this duration is selected. 
This interval is 180 points long and its R peak is located at 
the centre of the interval. This interval is selected as the 
signal of interest for matching the wavelet. 

The length of the wavelet is set to N=8, the average of 
the lengths of the wavelets used in previous works. The target 
function is the approximation error at the second scale (J=2). 
Following the steps discussed in the previous sub-section, the 
wavelet matched to this cardiac cycle is designed. The 
impulse response and frequency response of the 
corresponding low-pass filter, h, are shown in Fig. 1. 

4) Wavelet-Based R Peak Detection Algorithm 

The algorithm proposed in [1], [2] is the basis of R peak 
detection in this work. The algorithm is applied directly to the 
digitized signal without any filtering or pre-processing. This 
algorithm contains three main parts: 1) windowing the ECG 
signals to analyze 720 data points at a time, 2) computing the 

DWT of each interval at the scales 7} and 2' , and 3) 
analyzing the wavelet coefficients for R peak detection 

According to signal detection tteory, if the system 
response of the detection filter matohes the signal emtedded 
in noise, the signal-to-noise ratio (SNR) is maximized and a 
sharp and hic^ peak will be produced at the output showing 
the maximum correlation [12]. So using the matehed 
wavelets for deconposing the ECG signal generates sharp 
positive maxima- negative minima at the locations of R peaks 
at different scales as shown in Frg. 2. 

Concerning the power spectra of the ECG signal, rx)ise 
and artifact and the 3-db bandwidth of the equivalait wavelet 
anal5^s filter and in order to decrease the volume of 
mathematical conpitations for generating DWT coefficients 
at different scales, it seemed suffidait to analyze the wavelet 
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coefficients only at the scales 2^ and 2' . 

The R peak detection algorithm caii be divided into the 
following steps: 

Stepl: The DWT of the signal is calculated at scales 2^ 

and2' . 

Step 2: All the maxima in DWT larger than a threshold, 

winch is 1.7 times larger than the rms of DWT coeffiderts at 

each interval, are detected at scale 2 ^ . 

Step 3: All of the maximum-minimuin pairs in DWT larger 

than a threshold at scale 2^ are found. Similar to step 2, the 

threshold is updated after each detection. 

Step 4: In the neighboring of each maximum-minimum pair 

detected at scale, a similar pair is detected at scale 2' . 

Step 5: Isolated peals at scale 2' are eliminated. 

Step 6: Redundant peaks at scale 2' are eliminated. 

Step 7: The remaining pairs at scale 2' correspond to the R 
pedffi of the signal. The zero crossing points of these pairs are 
calculated and shifted 10 points to the left to oliain the 
locations of the R peals in the ECG signal. 

in. Results 

The proposed algorithm is evaluated on the standard and 
manually anrx)tated MIT-BIH Arrt^lhmia database. No pre- 
processing of filtering is ^jplied and the signals are directly 
analyzed using the proposed R peak detection algorithm 

To conpare the efficiency of the designed wavelet with 
the previously used wavelets, the second experiment is 
performed using quadratic spline wavelet \\irich had led to 
the best results among those of the previously proposed 
wavelet-based methods. In order to obtain a meaningful 
conparison among the different wavelets, the R peak 
detection algorithm remains mxhanged throuc^ tiiese 
experiments. 

Table 1 shows the results of the experiments for the 
MIT-BIH Arrhythmia data-base. 'TF' is the munber of false 
positives (incorrectly detected) and "FN" is the number of 
false negatives (missed) R peaks. The total number of failed 
detections (FD) is the sum of FP and FN. 

IV.CONCLUSION 

This work presented a new method for wavelet-based R 
peak detection by using matched wavelets. The previous 
works required analyzing the wavelet coeffidats in 4 
different scales vAale the proposed method needs analyzing 



ECG signal 









Fig. 1. Impulse response and frequency response of the low-pass filter which 
corresponds to the matched wavelet. 



Fig. 2. ECG signal and its DWT at scales 2 and 2 . The sharp maximum 
peaks at scales 2 and 2 correspond to R peaks. 

the coefficients of only 2 scales. 

The minimum execution time for the previous methods 
was about 108 seconds while for the proposed method it was 
reduced to 26.2 seconds. So, the speed of the algorithm is 
increased considerably. The results show a considerable 
increase in the accuracy of the R peak detection algorithm 
when using matched wavelet instead of Quadratic Spline 
wavelet, which had led to the best results among the previous 
wavelet-based algorithms. This is especially true for the 
number of miss detections, i.e., using matched wavelet 
decreases the probability of missing the R peaks notably. 
Future work will use further processing of the wavelet 
coefficients in order to increase the efficiency and reduce the 
effect of noise in the ECG signals. 
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TABLE I 

THE RESULTS OE THE R PEAK DETECTION FOR THE MIT/BIH DATABASE. RESULTS « CORRESPOND TO THE EXPERIMENT USING MATCHED WAVELET AND RESULTS t 

CORRESPOND TO THE EXPERIMNET USING QUADRATIC SPLINE WAVELET. 



Signal 


Total 
beats 


FP(#1) 


FP(*e) 


FN(#1) 


FN(*2) 


FD (#1) 


FD(ie) 


100 


2273 




















101 


1865 


2 


1 


1 


2 


3 


3 


102 


2187 








2 


3 


2 


3 


103 


2084 








4 





4 





104 


2230 


181 


137 


52 


57 


233 


194 


105 


2572 


24 


21 


32 


26 


56 


47 


106 


2027 





1 


101 


100 


101 


101 


107 


2137 


94 


651 


5 


486 


99 


1137 


108 


1763 


433 


525 


113 


132 


546 


657 


109 


2532 


4 


3 


7 


184 


11 


187 


111 


2124 


5 


8 


4 


6 


9 


14 


112 


2539 


2 











2 





113 


1795 








1 





1 





114 


1879 


6 


9 


7 


9 


13 


18 


115 


1953 




















116 


2412 


2 


13 


83 


48 


85 


61 


117 


1535 


1 


2 


1 


1 


2 


3 


118 


2275 


19 


4 


8 


3 


27 


7 


119 


1987 








4 


445 


4 


445 


121 


1863 


4 


2 


8 


2 


12 


4 


122 


2476 




















123 


1518 








3 


2 


3 


2 


124 


1619 


3 


33 


5 


4 


8 


37 


200 


2601 


108 


118 


174 


64 


282 


182 


201 


1963 


9 


9 


26 


66 


35 


75 


202 


2136 


1 





140 


162 


141 


162 


203 


2982 


101 


124 


168 


119 


269 


243 


205 


2656 


4 





32 


33 


36 


33 


207 


1862 


298 


317 


56 


90 


354 


407 


208 


2956 


23 


19 


422 


642 


445 


661 


209 


3004 


11 


6 


1 


18 


12 


24 


210 


2647 


30 


14 


38 


77 


68 


91 


212 


2748 


1 


2 





3 


1 


5 


213 


3251 


5 


1 


12 


175 


17 


176 


214 


2208 


18 


8 


16 


26 


34 


34 


215 


2154 


11 


9 


8 


16 


19 


25 


217 


2208 


11 


14 


16 


42 


27 


56 


219 


2154 


4 


6 


2 


26 


6 


32 


220 


2048 








23 


30 


23 


30 


221 


2427 


3 


1 


14 


24 


17 


25 


222 


2484 


12 


16 


44 


58 


56 


74 


223 


2605 








9 


11 


9 


11 


228 


2053 


363 


168 


40 


21 


403 


189 


230 


2256 


19 





2 





21 





231 


1886 








1 


1 


1 


1 


232 


1780 


793 


588 


1 


5 


794 


593 


233 


3079 


43 





57 


71 


100 


71 


234 


2753 











2 





2 


total 


116137 


2648 


2830 


1743 


3292 


4391 


6122 
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